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We argue that hadron production in coherent diffraction of proton on a heavy nucleus 
provides a very sensitive probe of the low- a; QCD dynamics. This process probes the BFKL 
dynamics in proton and the non-linear gluon evolution in nucleus. We calculate the diffractive 
hadron production cross sections in the RH1C and LHC kinematic regions. To study the 
nuclear effects we introduce the diffractive nuclear modification factor. We show that unlike 
the nuclear modification factor for inclusive hadron production that has a very interesting 
dynamics at RH1C but is expected to be almost completely saturated at the LHC, the nuclear 
modification factor for diffractive production exhibits a non-trivial behavior both at RH1C 
and LHC. 

I. INTRODUCTION 

Over the last decade we have witnessed a remarkable success of models based on gluon saturation 
in description of low- a; data at HERA and RHIC. This allowed quantification of several key features 
of the low-x dynamics of QCD. Still, there are many open questions such as the size of the NLO 
corrections to the BK equation, demarcation of boundary between the kinematic regions of gluon 
saturation and the collinear factorization, etc. These problems can be addressed by probing the 
nuclear structure at even smaller x and/or by using a different set of measurements. In this paper 
we argue that diffractive hadron production in pA collisions is a measurement that can provide a 
new handle on the low-x nuclear dynamics. Our study is motivated by the possibility to investigate 
the diffractive processes using the data on deuteron-gold (D-Au) collisions collected at RHIC. 

A detailed theoretical analysis of coherent diffractive gluon production in onium-heavy nucleus 
(qqA) collisions in the framework of the Color Glass Condensate [U [2j (3J IH |5j [6] was performed 
in our previous publications |5]. There we argued that this process is sensitive to the low-x 
dynamics both in onium and nucleus. It has been argued in [9], [TU1 [TTj that it is phenomenologically 
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reasonable to approximate the proton light-cone wave-function (away from fragmentation regions) 
by a system of color dipoles. Additionally, we will demonstrate in Sec. [IT|that the qqqA propagator 
in the quasi-classical approximation takes exactly the same form as the qqA one if the emitted 
gluon transverse momentum is hard kx 3> Q s . In this approximation we can directly adopt the 
results of our theoretical analysis in Ref. |8] . The corresponding phenomenological approach is 



developed in Sec. Ill Similar model has been used in Ref. |12j for description of diffraction in pA 
collisions. 

There are several parameters that govern behavior of diffractive gluon production in pA colli- 
sions. These are: gluon transverse momentum kx and rapidity y, nucleus atomic number A and 
transverse distance between the valence quarks rx in proton. The main observation of 18] is that 
dependence of the diffractive hadron spectrum on these parameters in various kinematic regions is 
quite different. This provides a convenient handle on the behavior of the low-x gluon densities in 
the three most interesting kinematic regions: (i) gluon saturation region kx < Q s , (h) geometric 
scaling region kx < Qgeom and (hi) hard perturbative QCD region kx > Qgeom- 

The model that we use in this paper is based on analysis of diffractive hadron production in all 



available kinematic regions. Eq. (21 ) holds in the logarithmic approximation in all those regions and 
is therefore a convenient interpolation formula which we use to calculate the differential inclusive 



cross section (15). The transverse vector I(r',k, y) encodes information about the gluon density 
in the nucleus. It is related to an integral of the forward elastic gluon dipole scattering amplitude 



Na(t' ,h,y) over all intermediate dipole size, see (16) and (17). This amplitude is parameterized 
according to the KKT model [13 . On the other hand, the dipole density n(r, r', Y — y) encodes the 
gluon density in the proton, which is assumed to be dilute. Since the dipole density is a solution 
to the BFKL equation, we model it by the LO BFKL amplitude in the diffusion approximation. 
To compare the low-x dynamics in pA collisions to that in pp ones, it is convenient to intro- 



duce the diffractive nuclear modification factor R^, see (31). We evaluate the diffractive gluon 
production in pp collisions as a limit A — > 1 of that in pA ones. Theoretical expectations for R pA 



are detailed in Sec. Ill B and Sec. Ill C The results of our numerical calculations performed using 



the KKT model [13] are presented in Sec. IV[ We observe, that R^ behavior is quite different 



from that of the nuclear modification factor R? , for inclusive hadron production. In the RHIC 
kinematic region, at moderately large kx there is a significant enhancement of particle production 
in pA collisions, see Fig. [3} This happens due to the fact that the diffractive cross section at large 
kx is proportional to the higher twist contribution that is enhanced in pA collisions by an addi- 
tional factor of A 1 ^. This enhancement gets increasingly compensated at forward rapidities by a 
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suppression stemming from two sources: (i) gluon saturation in the nucleus; (ii) shrinking of phase 
space available for the BFKL evolution in proton jj] . The latter feature of the diffractive hadron 
production is apparent in (34) and is illustrated in Fig. [5] where we compare R^g for two different 



diffusion coefficients (switching the BFKL evolution on and off). In Fig.|4]and Fig. [(3]we show R p ^ s 



Miff 

at LHC. Rfoff exhibits rather strong dependence on rapidity. In contrast, Rf nc i is not expected 
to change a lot at LHC |15j. This implies, that by comparing inclusive and diffractive hadron 
production in the wide kinematic region of RHIC and LHC one will be able to infer much of useful 
information about the higher twist contributions. Since different models of low-x dynamics predict 
different dependence of higher twists on atomic number A and energy /rapidity, measurements of 
diffractive hadron production will be instrumental in determining the valid physical mechanism for 
hadron production at high energies. 

II. A MODEL FOR DIFFRACTIVE GLUON PRODUCTION IN PA COLLISIONS 
A. Diffractive gluon production in qqqA collisions 

Coherent diffraction of a proton on a nucleus is a process p + A — > X + A characterized by a 
large rapidity gap between the diffractive system X and the intact nucleus A. A fraction of the 
coherent diffractive events increases with the collision energy and is expected to reach its limiting 
value of a half at asymptotically high energies. In the mean- field approximation a s <C 1 and A>1, 
the incoherent diffractive processes such as p + A — > X + A* where A* is a diffractive system of 
color-neutral nuclear debris, are parametrically suppressed. Therefore, in the present paper we 
consider only the coherent diffraction 1 . 

Coherent diffraction is possible only if the coherence length l c of the emitted gluon with mo- 
mentum k is larger than the nucleus size Ra (in the nucleus rest frame): 

Ic = ^^Ra, (1) 

where + indicates the light-cone direction of the incoming proton. The invariant mass of the 
produced system is given by M 2 = k 2 /x, where x = /c+/p+ and p is the proton momentum. 
Substituting these equations in ([I]) yields the following condition on the mass of the diffractive 
system: 

M 2 « ^+ = — , (2) 
Ra Ra^p 



1 The incoherent diffraction may be phenomenologically important at RHIC and LHC energies [121 116j . 
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where y/s is the center-of-mass energy of the proton-nucleon collision and m p is proton mass. 

A realistic model for diffractive gluon production in pA collisions was discussed by Kovchegov in 
|14| . He considered, in the quasi-classical approximation, emission of a gluon by a color-neutral qqq 
system of valence quarks with subsequent elastic interaction with a heavy nucleus. The resulting 
expressions for the propagators of the qqq and qqqG systems in the nucleus can be written in the 
following form [J2] : 

Uij = vf (V M(Zl) - e -( 2 /9)(xi2+Xl3+X23)^ L-M( Z2 ) _ e -(2/9)( X i2+Xi3+X23)^ v . (3) 

where 

vl = (-1, ~), vl = (1, -i=), vl = (0, -|) , (4 ) 



the 2x2 matrix M(z) is given by 
M(z) = 



gC3 + ^(C2 + Cl) + ^(X23 + X13) - |X12 175 (~^ 2 + Cl + X23 ~ X13) 

J^(-C2 + Cl + X23 - X13) |C3 + l(C2 + Cl) - m(X23 + X13) + §X12 



(5) 

and the scattering amplitudes of various dipoles on a nucleon read 2 

Cf = - x i) 2 Q 2 sO Xij = j^(xj - KjfQlo > ( 6 ) 

where xi, X2, X3 are the valence quarks transverse coordinates, zi and Z2 are the gluon transverse 
coordinates in the amplitude and in the complex conjugated one respectively, see Fig. [TJ 

The cross section for the diffractive gluon production in the quasi-classical approximation reads 

& it w 

Assume that the distances between the valence quarks are approximately the same: Xij ~ X = 
|-RpQ 2 Q, where R p is the proton radius. Then, matrix M from ^ reduces to 

M(z)=^ C3 + A (C2 + Cl) + ' X iJl^ + W "| (8) 

V j^<-<2 + 0) |C3 + J(C2 + Ci) + |x/ 

In general, the propagators Ilj,- are rather complicated objects. However, in the perturbative regime 
depicted in Fig. [T^a) they can be reduced to a simple sum of the corresponding qqG propagators as 



2 In this section only we adopted a shorthand notation where the saturation scale is understood to include the 
logarithmic dependence on the dipole size. 
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FIG. 1: Possible configurations formed by three valence quarks at positions xi, X2, X3 and a gluon at 
position Zj. We assumed that X12 ~ X23 ~ 2:31- In (a) sizes of the daughter dipoles are much bigger than 
sizes of the parent dipoles. For small enough X12 this corresponds to the pQCD regime. In (b) one of the 
daughter dipoles is much smaller than the rest of dipoles corresponding to the high density regime (see [T7] 
for more details). 



we are going to demonstrate now. The leading logarithmic contribution in the perturbative regime 
stems from the configuration shown in Fig. [ija). In this case Q S> x an d we have Ci ~ C2 ~ C3 = C- 
Eq. ^ becomes 



M(z) 



Using ^ and @ we derive the propagator 




(9) 



_ 2 

e 3 



\ \ 1 e -(M 



e 3 



(10) 



The cross section ^ reads in this case 



d 2 k T dy ~ (2vr) 2 
3 



j d 2 bd 2 Zl d 2 z 2 e - k -( Z i- Z2 ) 



1 



E 



Zi - Xi 



Zi - x 7 - 



2 f—r \|zi - Xj 



z l — x 



■ |2 



Z2 ~ Xj 
|z 2 - Xjl 



z 2 - x,- 



Z2 - X 



■ |2 



n 
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(ii) 



In the 't Hooft's limit, Eq. (11) can be related to the cross section for diffractive gluon production 



in quarkonium-nucleus collisions. To this end we introduce an effective color dipole with a quark 
and antiquark being at points xi and X2 respectively. Then, in the same approximation as in 
Fig. [ija), we obtain 

do qq ~ A ^ o. s Cf 1 f ,2, ,2 ,2 ( z i — x i z i ~ x 2 \ ( z 2 - xi z 2 - x 2 
d 2 kr dy ir 2 (2tt) 2 



d 2 bd 2 Zl d 2 z 2 



zi - Xi 



zi - x 2 



|z 2 -Xi| 2 |Z2-X 2 | 2 



x e -ik-(zi-z 2 ) ( e -\{x\-*l) 2 Q 2 s0 



e s 



i(xi-x 2 ) 2 Q2 Q \ / _l (il _ Z2 )2^ 



;-(xi-x 2 ) 2 



Hence 



da qqqA 2 da 9 ^ 3 da qqA 



d 2 kxdy Cf d 2 k^dy 2 d 2 kxdy 



(13) 



6 



Comparing arguments of exponents in (10) and in (12) we identify Q 2 = \Q 2 S $ as an effective 
saturation scale and 

R 2 = (X! - x 2 ) 2 = 2 • | ( Xl - x 2 ) 2 = \ ■ 3 R 2 = (2 i? p ) 2 (14) 



as the square of the dipole size. Expression (13) motivates a model that we adopt in this paper. 



We assume that the pA cross section can be approximated by qqA one with the dipole size given 



by (14). This model correctly reproduces the pQCD limit. It also satisfies the unitarity bound, 



which is achieved in the saturation regime depicted in Fig. [T]^b). 

B. Gluon production in quarkonium heavy nucleus collisions 

Now, as we set up a model for the diffractive gluon production in pA collisions in terms of the 
diffr active gluon production in qqA collisions, we would like to review the main results that we 
derived for the latter case in our previous publications [3 |8] . The cross section for the diffractive 
gluon production with transverse momentum kj> at rapidity y is given by 

where n p (r,r',Y — y) is the dipole density in the projectile proton. It has the meaning of the 
number of dipoles of size r' at rapidity Y — y generated by evolution from the original dipole r 



having rapidity Y [18]. It satisfies the BFKL equation |19| [20] with the initial condition (18). The 
two-dimensional vector function I(r',k, y) is defined as follows: 

I(r', k, y) = -e- iW zV k Q(r', k, y) + iV k Q*(r', k, y) , (16) 

where 

Q(r',k,y) = - / d 2 we^\ 
J w l 

x [N A {r', b, y) - N A (w - r', b, y) - N A (w, b, y) + N A {w - r' , b, y)N A (w, b, y)] . (17) 

The vector function I(r',k, y) incorporates information about two physical processes: (i) gluon 
emission off the daughter dipole r' produced in the course of the BFKL evolution and (ii) low- 
x gluon evolution in the nucleus through N A (r,h,y), which is the dipole-nucleus forward elastic 



3 In the following we are going to discuss only onium-nucleus scattering. Therefore we will omit the tildes to simplify 
notations. 
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scattering amplitude satisfying the BK equation |21[ [22j. In the quasi-classical approximation the 
dipole density reads 



n p (r, r', 0) = 5(r — r') 



(18) 



while the scattering amplitude is given by the Glauber-Mueller formula [23] (now we explicitly 
write down the logarithm in the exponent) 



JV A (r,b,0) = 1 - e -t r2Q 'oi°^ 



(19) 



where Q s o is the saturation scale at rapidity y = 0. In the case of dipole-proton scattering we 



expand (19) and get 



iV p (r,b,0) = ^r 2 A 2 ln-i-. (20) 
In all limiting cases we can write [8] 

|I(r', k, y)\ 2 « C N \{k-% b, y) [1 - N A (r', b, y)] 2 sin 2 , (21) 

where C is a constant of order unity (its precise value which can be found in |8j is of little importance 



here). Let us emphasize, that (21 ) holds asymptotically in all kinematic regions. Due to the initial 



condition (18), the cross section in the quasi-classical approximation is merely proportional to 



|I(r', k, y)\ 2 . Accordingly, employing (21) we obtain 

da pA {R,k T ,0) 
(Pkxdy 



C 4 ^^N 2 Ak T \h 1 0)[l-N A (RM0)] 2 sin 2 (^^) , (22) 



At larger rapidities we integrate over r' in (15 ) using (21 ). In the case of hard gluons we get 



dgP A (R,k T ,y) 
d 2 kxdy 



CXsCf 



S A Ni(k-\b,y) 



min^i? 2 } 



,2y/2a s (Y-y)\ln(Rk T )\ 



, k T > Q s , (23) 



775/2 ' ' (2a s (Y-y)\ln(Rk T )\) 1 / 4 

where a s = a s N c /TT. The cross section for the soft gluon production by a large dipole reads 



da pA (R,k T ,y) = a s C F Sa {2a s {Y - y)) l l 4 2^2a s (Y-y)HRQ. 

d*k T dy 87T5/2 Q2 ln 3/4 (i?Qs) 

while in the case of soft gluon emission by a small onium 



do-P A (R,k T ,y) a s C F 



d 2 krdy 



4vr 5/2 



S A R 2 



1 



2a s (Y-y)ln 1 i 



1/4 



2^/2a s (Y-y) In 



TT > 7T ' (24) 
k T Qs 



i?<i-< (25) 

Qs K T 



In all the reviewed cases (23)-(25) gluon multiplicity arises from the cut Pomeron that is hooked 



up to the incoming proton. 
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C. Forward dipole nucleus scattering amplitude 

The last required ingredient is the forward elastic scattering amplitude Na(t, b, y). It can be 
evaluated in various kinematic regions. In the double logarithmic approximation (DLA) 

In 1 / 4 ( -4- 



r<l/Q s0 , In— ^->a s y. (26) 

This limit coincides with the small x and small r limit of the DGLAP equation. It obviously 
breaks the geometric scaling. Consequently, the DLA holds in the transition region between the 
gluon saturation and the hard perturbative QCD characterized by a hard scale kn , i.e. when 
Qgeom < kx < kjj, where Q ge0 m is the scale at which the geometric scaling breaks down. It reads 
in the DLA 

Ql 

Qgeora ~ 7^ • (27) 
VsO 

The saturation scale is given by 

Q s ps A 1 '* A 2 e XY , (28) 
where A ps 2a s in the DLA. The hard scale kn can be related to the invariant mass of the 



diffractively produced system as discussed in Sec. II A in detail. 

As we approach the saturation region by decreasing kx at fixed rapidity we arrive at the diffusion 
approximation 

"^-Iti/toKt)^'"^ — 1,,2 (^)' (29) 

where the BFKL Pomeron intercept is ap — 1 = 4a s ln2. We observe that the amplitude geo- 



metrically scales modulo small diffusive corrections. The diffusion approximation (29) holds in the 
kinematic region Q s < k? < Qgeom- 

Finally, when k? < Q s solution to the BK equation deeply in the saturation region implies [T7] 

N A (r,b,y) = l-S e-T 2 /* = l-Soe~> 2 ^\ r > -L . (30) 
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III. NUCLEAR EFFECTS IN DIFFRACTIVE GLUON PRODUCTION 
A. Nuclear modification factor 

A convenient way to study the nuclear dependence of particle production is to consider the 
nuclear modification factor defined as follows 

d<J diff( fc T,J/) 



d 2 kj-dy 



If the production process is completely incoherent then R P ^{kx, y) = 1. In the case of inclusive 
gluon production, the nuclear modification factor i?? , was discussed in detail in |24j . It has been 
demonstrated that in the extended geometric scaling region Q s (y) kx < Qgeom, the nuclear 
modification factor is suppressed as Rf^i ~ j4 _1//6 , while in the saturation region kr < Qs(y) 
the suppression is 

^fnd ~ A' 1 / 3 . The amount of suppression is closely related to the value of the 
anomalous dimension 7 in a given kinematic region. At rapidity y ~ at RHIC Rf nc \ exhibits slight 
enhancement (Cronin effect), which serves as indicator that the low-x evolution in that process 
does not play an important role. We are going to argue below that the behavior of R%ig is quite 
different from that of inclusive one which makes it a convenient tool for study of the low- 2; gluon 
dynamics. 

In the previous section we addressed in detail the diffractive gluon production in pA collisions. 
In order to evaluate the R^ iS we need to normalize it by that in pp collisions. The latter is 



obtained by replacing the forward elastic dipole-nucleus scattering amplitude given by (26), with 
the corresponding forward elastic dipole-proton scattering amplitude 



In (26) we replaced Q s q by A and set A = 1. Since we assume that the gluon saturation effects are 



negligible in the proton, the cross section for the diffractive gluon production in pp collisions in 



the case of large characteristic proton size is obtained from (23) by setting A = 1 with the result 



d 2 k T dy 7T 5 / 2 14 j p P y t (2a s (Y - y)\ln(Rk T )\) 1/4 

(33) 



Similarly, we get in the quasi-classical approximation using (22) 



^ Wt ,0, , c <*|^ k ,(M H^'id-JbtfUr)), (34) 



d 2 k T dy ir 2 k 2 64 4 \A 
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where we averaged over the directions of the dipole R according to 

i f dO sin 2 ( 1 k T Rcos6) = 1 (1- J (Rk T )) • (35) 
^ Jo V 2 / 2 

Gluon saturation effects in proton may be important at backward rapidities at the LHC. Taking 
them into account constitutes a difficult and not yet solved problem. Fortunately, effects associated 
with the gluon saturation in proton are not expected to significantly alter the nuclear dependence 



of our results since they are likely to cancel between the numerator and denominator of (31) 



B. Quasi-classical approximation 

The nuclear modification factor in the quasi-classical approximation and at high transverse 



momenta is derived by substitution of (22) and (34) into (31) and deduce 



<(k T , 0) = AW (l - \a 1/3 ^ In ^) e-i* 2 ^ , kT » Q sQ , (36) 



where we take into account that Sa = A 2 ^ 3 S P and Q 2 = A 1 / 3 A 2 . According to (36) at very large 
and fixed A the nuclear modification factor approaches a constant 

< ff (^T,0)-^/V^ 1/3l ^ 1/3 , k T ^oo. (37) 



Eq. (36) implies that R^g (&T 5 0) approaches unity from below as kT — ► oo. In contrast to 



R P J[R(k,Ti 0), the nuclear modification factor for inclusive gluon production receives a positive power 
correction that is a source of the Cronin enhancement observed in inclusive gluon production in 
pA collisions. 

In the saturation region we derive 

^ fe0) = iiWX e " W - k T «Q s ,. (38) 
That is, the nuclear modification factor vanishes at small momenta as k T . Actually, if we neglect the 



slow logarithmic dependence of the initial saturation scale Q s q on r in (19) the integral appearing 



in (17) can be taken analytically. The corresponding result can be found in p3]. In Fig. [2] we 
use this analytical result to plot the nuclear modification factor R p ^ as a function of transverse 
momentum kT- 

We observe that, unlike in the inclusive gluon production case, the size of the incoming projec- 
tile plays a very important role in the diffractive production. What is important is the relationship 
between the quarkonium size R and the inverse saturation scale 1 /Q s . In the quasi-classical approx- 
imation, that is neglecting the low-x evolution, the diffractive gluon production is exponentially 
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FIG. 2: Nuclear modification factor R^ is as a function of transverse momentum kx in the quasi-classical 
approximation. A is a non-perturbative momentum scale. 

suppressed for heavy nuclei if R > 1/Q S as compared to light nuclei. If R < 1/Q S suppression 
gives way to enhancement at high transverse momenta. Both effects come about as the result of 
the coherent scattering of proton off nucleus. 



C. Low- a; evolution: hard gluons 



1. Double logarithmic approximation 



The low- a; evolution has a dramatic effect on the diffractive gluon production. We would like to 
start our analyses with the case of moderately large transverse momentum such that the geometric 



scaling is broken, but interaction is still coherent. Substituting (23) and (33) into (31) we derive 
that in general 



R p d f s (R,k T ,y) 



1 N\{k T \h,y) 



k T ^>Q s . 



(39) 



A^Ni(k T \b,y) 

In the double logarithmic approximation, the BFKL equation coincides with the DGLAP equation. 
Therefore, in this region we can observe crossover from the coherent small- x dynamics to incoherent 



hard perturbative QCD. Using (|26|) and (|32|) in (|39|) we derive 
S A 



R p d tf(k T ,y) 



A Snr> 



lni^^i 



6-0 



Qso Q 



1 + 



2a s y . Qso 



4 v / 2aly 



In 

Qso 



In 



kx S> Qgcom • 

(40) 
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Introducing a new variable 



(41) 



we reduce ( 40 ) to 



exp < -4W 2a s y In 



Qso i-C 2 



a y 1 + c 2 



(42) 



The DLA approximation is valid when kx 3> Qs(y) > Qso- I n this case (42) becomes 



R P A(kr,y) » AV3 1-^2 



In ^ 
m A 



k T > Q P 



(43) 



The remarkable feature of this result is enhancement of the nuclear modification factor by A 1 / 3 . 



Unlike the quasi-classical case (37), this enhancement is not overrun at large A by a small exponen- 



tial factor. The reason is that in course of low-x BFKL evolution dipoles with small size r < 1/Q S 
are produced and these dominate the cross section. Let us also mention that an enhancement 



similar to (43) has already been discussed in context of the J/i/j production off the nuclear targets 
|27j as well as in the breakdown of the collinear factorization of the fragmentation functions |28j. 



It is important to emphasize that the result (43) holds only as long as the coherence length 



lr 



I — - is much larger than the nuclear size. Since in the center-of-mass frame kinematics 



2M N x 



at large enough transverse momentum fey and fixed rapidity y and energy s the 



coherence is lost and the nuclear modification factor approaches unity. Therefore, the region where 
R dfs ~ A±/3 scaling gives way to -R^iff ~ 1 is the transition region between the semi-hard nuclear 
fields and the hard perturbative QCD. Needless to say that identification of this region is crucial 
for understanding the interplay between the dense and dilute high energy QCD regimes. 



2. Extended geometric scaling region 

Next, we would like to analyze the extended geometric scaling region Q s {y) < kx < Qgeom- 
Here the evolution is still linear and is well approximated by the leading twist approximation. 
However, the anomalous dimension of the gluon distribution significantly departs from unity and 
approaches the value it has at the critical line kx = Qs{y)- It is therefore appropriate to use the 



leading logarithmic approximation for the function iV^r, b, y). Substituting (29) and (32) in (39) 
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we derive 

vA n \ 4 &h?(%e)y/2ad f / 21n 2 (%° 



7C(3)A2 -r^-r ^ * V ~„~ A 14C(3) ^ 

geom • 

(44) 

This equation clearly demonstrates that the ^-dependence of the nuclear modification factor arises 
only through the slow- varying logarithmic factors. As far as the rapidity dependence is concerned, 
we can estimate it at the scale kx = Qgeom(y)- Since Nj^{r, b, y) is constant on the critical line we 
derive 

R p d f s (Q gcoia (y),y)~A 1 ^e-^y. (45) 

That is, the nuclear modification factor is getting progressively suppressed in the forward direction. 
This is much stronger suppression than in inclusive gluon production. Approximately we can write 

Rm(kT,y)-A 1 / 3 (^ cl (k T ,y)) 2 , Q s <k T <Q gcom . (46) 



Eq. (46) clearly exhibits the higher twist nature of the diffractive gluon production. The peculiar 
properties of diffractive cross section due to the higher twist contributions in nuclear and hadronic 
DIS have been discussed in [25j ES] . 

3. Saturation region 



In the saturation region we utilize (33) and one of the (24) or (25) in (31) and arrive at a rather 
involved expression. Keeping only the parametric dependence and omitting the logarithmic factors 
we obtain 

R P A(R, k T , y) ~ ^jSq-^ 2 ^ 3 ^ 1 ^ e- 2 ^MV-,)m(*M e -4vW^ , kT , I < Q s 

(47) 

There is a very strong suppression of diffractive gluon production in the saturation region in the 
case of low-x evolution. This suppression however is still milder than in the quasi-classical case 
(|38[). On the critical line kx = Qs(y) we get for forward rapidities (Y — y<^y) and central collisions 



(employing (28)) 

R p A(Qs(y),y) - e-^y , (48) 

which implies a strong suppression in the forward direction. 
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IV. NUMERICAL CALCULATIONS 

All the features that we discussed in the previous section can be visualized using a simple 
model for the forward elastic dipole-nucleus scattering amplitude iVyi(r, b, y). We parameterize it 
as follows [13] 

iV A (r,b,y) = l-exp|-^(r 2 g2)7(^)| . (49) 

The anomalous dimension is parameterized in such a way that it satisfies the analytically well- 
known limits of (i) r — > 0, y fixed and (ii) y — > oo, r fixed: 

i (i + , 4jsa^ r —— > \ „ > yo 



where 



i y <yo, 



In [l/(r 2 Q 2 )] 

( A / 2 )(y - Vo) 

In the double logarithmic approximation we can replace r 2 ~ 1/(4A;|>). The gluon saturation scale 
is given by 

«-« = a ^ 1/v "(m) X < < 52 > 

where parameters A = 0.6 GeV and A = 0.3 are fixed by DIS data [30 . The initial saturation scale 



used in (51 ) is defined by Q sQ = Q s (yo) with y$ the value of rapidity at which the small-x quantum 



evolution effects set in. Fit to the RHIC data yields yo = 0.5 |13j . 



Numerical calculations of the cross section (15) are performed after substitution of (21) with 



(49) and the following formula for the dipole density in diffusion approximation (cp. (29)): 



Mry,Y-v) = _ g ^-" y -^-"^^ ■ (53) 

Parameter d is equal to unity in the LO BFKL. To obtain the hadron diffr active cross section we 
convoluted the obtained result with the LO pion fragmentation function given in (29]. Diffractive 
gluon production in pp collisions, which is required as a baseline for the calculation of the nuclear 



modification factor (31), is obtained by setting A = 1 in the formula for the corresponding cross 
section in pA collisions. 

The results of numerical calculations are exhibited in Fig. [3]-Fig. [6] In Fig. [3] one can see that at 
RHIC R pA ~ 2 - 2.5 at y ~ and k T > 2 GeV. This enhancement is a signature of a leading power 
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FIG. 3: Nuclear modification factor for diffractive pion production in pA collisions at RHIC as a function 
of transverse momentum for two characteristic sizes of proton (a) 0.2 fm and (b) 1 fm. The effects of finite 
coherence length are neglected. 



r=l GeV" 1 , Vs~ =5500 GeV, d=l 




(a) 



(b) 



FIG. 4: Nuclear modification factor for diffractive pion production in pA collisions at LHC as a function of 
transverse momentum for two characteristic sizes of proton (a) 0.2 fm and (b) 1 fm. The effects of finite 
coherence length are neglected. 



correction, see (43). As rapidity increases there are two important effects, which take place in 
the proton and nucleus wave functions: (i) spectrum of intermediate dipoles in a projectile proton 
shrinks as the rapidity interval available for the low-x evolution in proton becomes narrower, (ii) 
as y increases, x of gluon decreases causing stronger gluon saturation effect in the nucleus. Both 
effects lead to suppression of the nuclear modification factor. Gluon saturation in proton leads 



to the suppression law (45). Of course, the effect of diffusion in a proton is more pronounced 



for a proton with larger characteristic size, since in absence of the evolution effects (i.e. in the 



quasi-classical approximation) the cross section would be exponentially suppressed, see (38). 



We further investigated the effect of diffusion by introducing the parameter d in (53). As has 
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FIG. 5: Effect of diffusion in the dipole sizes on the diffractive pion production at RHIC for two characteristic 
sizes of proton (a) 0.2 fm and (b) 1 fm. Upper line of the same type corresponds to d = 1, the lower one 
- d — 0.1. Lines of different types correspond to different rapidities (notations are the same as in Fig. [3]). 
The effects of finite coherence length are neglected. 
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FIG. 6: Effect of diffusion in the dipole sizes on the diffractive pion production at LHC for two characteristic 
sizes of proton (a) 0.2 fm and (b) 1 fm. Upper line of the same type corresponds to d — 1, the lower one 
- d = 0.1. Lines of different types correspond to different rapidities (notations are the same as in Fig. [1]). 
The effects of finite coherence length are neglected. 



been repeatedly pointed out in this paper, it is the BFKL diffusion that makes the diffractive gluon 
production possible by generating intermediate dipoles of small size. Gluon saturation effects in 
proton may tame the BFKL diffusion (31] leading to smaller effective diffusion coefficient. This 
effect is taken into account in Fig. [5] for RHIC and in Fig. [6] for LHC. The shadow region in all 
figures demonstrates the difference in the nuclear modification factor between the cases of d = 1 
and d = 0.1. Switching off the diffusion severely impacts the nuclear modification factor at low 
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energies/rapidities and for larger distances between the valence quarks in proton. 

Unlike the nuclear modification factor for inclusive hadron production, which decreases as a 
function of rapidity and centrality at RHIC and reaches almost maximal possible suppression so 
that no significant additional suppression is expected at LHC, the diffractive hadron production 
shows a very interesting behavior even at LHC. This makes this process suitable for exploration of 
different kinematic regions at the high energy frontier. We believe that it will be instrumental in 
unraveling the structure and dynamics of strong gluon fields. 
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